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We present a careful and thorough microscopic derivation of the Kondo Hamiltonian for single-molecule 
magnets (SMMs) transistors. When the molecule is strongly coupled to metallic leads, we show that by applying 
a transverse magnetic field it is possible to topologically induce or quench the Kondo effect in the conductance 
of a SMM with either an integer or a half-integer spin S > 1/2. This topological Kondo effect is due to the 
Berry phase interference between multiple quantum tunneling paths of the spin. We calculate the renormalized 
Berry phase oscillations of the two Kondo peaks as a function of the transverse magnetic field by means of the 
poor man's scaling. In particular, we show that the Kondo exchange interaction between itinerant electrons in 
the leads and the SMM pseudo spin 1/2 depends crucially on the SMM spin selection rules for the addition and 
subtraction of an electron and can range from antiferromagnetic to ferromagnetic. We illustrate our findings with 
the SMM Ni4, which we propose as a possible candidate for the experimental observation of the conductance 
oscillations. 
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I. INTRODUCTION 

Single-molecule magnets (SMMs), such as Mni2 (see Refs. 
[312) and Feg (see Refs. 3^), have become the focus of intense 
research since experiments on bulk samples demonstrated the 
quantum tunneling of a single magnetic moment on a macro- 
scopic scale. These molecules are characterized by a large 
total spin, a large magnetic anisotropy barrier, and anisotropy 
terms which allow the spin to tunnel through the barrier. Elec- 
tronic transport through SMMs offers several unique features 
with potentially large impact on applications such as high- 
density magnetic storage as well as quantum computing. 5 Re- 
cent experiments have pointed out the importance of the in- 
terference between spin tunneling paths in molecules. For in- 
stance, measurements of the magnetization in bulk Feg have 
observed oscillations in the tunnel splitting A m? _ m between 
states S z = m and — m as a function of a transverse mag- 
netic field at temperatures between 0.05 K and 0.7 K (see Ref. 
0). This effect can be explained by the interference between 
Berry phases associated to spin tunneling paths of opposite 
windings. 7,8 Theoretically, a coherent spin-state path integral 
formulation is used to account for the coherence of the virtual 
states over which the spin tunnels, although the initial and fi- 
nal spin states do not retain their coherence. 

A new approach to the study of SMMs opened up recently 
with the first observation of quantized electronic transport 
through an isolated Mni2 molecule. 9 One expects a rich inter- 
play between quantum tunneling, phase coherence, and elec- 
tronic correlations in the transport properties of SMMs. In 
fact, it has been argued that the Kondo effect would only be 
observable for SMMs with half-integer spin 10,11 and therefore 
absent for SMMs such as Mni2, Feg, and N14, where the spin 
is integer. Later, two of us showed that this is not the case: 12 
Remarkably, a transverse magnetic field H± can be tuned to 
topologically quench the two lowest levels of a full-integer 
spin SMM, making them degenerate. In fact, the same Berry- 
phase interference also affects transport for SMMs with half- 



integer spin: In that case, sweeping H± will lead not to one 
but a series of Kondo resonances. In the case of SMMs, as we 
show below, the Berry phase oscillations of the tunnel splitting 
^m ,m' of the uncharged single-molecule magnet (q = 0) 
for both full- and half-integer spins leads to oscillation of the 
Kondo effect as a function of H± . This means that the Kondo 
effect is observable at zero bias for all values of the magnetic 
field #_L, such that A mo ^ o (ifj_,o) = 0. 

It is interesting to note that at a finite bias the Kondo effect 
in a quantum dot in the presence of a magnetic field can be re- 
stored by tuning the bias to eV = ±gfiBH_\_ (see Ref.[l3|). For 
SMMs, however, the interference between the Berry phases of 
the molecule total spin makes the distance between the split 
Kondo peaks, which is equal to eV = ±A m0jm / , oscillate 
as a function of H±. A necessary condition for observing 
these oscillations is a large enough tunnel splitting. Recently 
a new SMM based on tetranuclear nickel clusters N14 with a 
S = 4 ground state has been synthesized. 14 Our motivation 
for studying this particular nanomagnet stems partly from its 
high symmetry (S4) but, more importantly, from the large tun- 
nel splittings: A mo5 _ mo ~ 0.01 K or larger (depending on the 
transverse field H±) between the \m) = |4) and \m') — \ — 4) 
ground states. 

Recently, some authors have argued that the Kondo ef- 
fect is absent at the diabolic points of the Berry-phase 
interference. 11 This conclusion came from considering an ad 
hoc Kondo HamiltonianiSJJ. which was not derived micro- 
scopically. Moreover, in recent analysis of the sequential tun- 
neling regime, an ad hoc exchange Hamiltonian mixed with an 
Anderson-type Hamiltonian was also used without a proper 
microscopic derivation. 15 In this paper, we provide a care- 
ful microscopic derivation of the Kondo Hamiltonian suit- 
able for full- and half-integer spin single-molecule magnets 
by means of a Schrieffer- Wolff transformation. 16 By using 
the exact eigenstates of the positively and negatively charged 
single-molecule magnet, it is sufficient to apply the Schrieffer- 
Wolff transformation to second order in the tunneling matrix 
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element. The resulting Kondo exchange parameters exhibit 
the interference between the second-order transition paths go- 
ing over the two virtual charged states. We show that this 
very same interference phenomenon is also responsible for 
the Berry-phase blockade of the current through the single- 
molecule magnet in the cotunneling regime, which extends 
our previous results obtained for the Berry-phase blockade in 
the sequential tunneling regime.— 

Our derivation of the Kondo Hamiltonian reveals an im- 
portant detail: The Anderson-type Hamiltonian of the single- 
molecule magnet can be mapped onto a spin- 1/2 Kondo 
Hamiltonian in two different ways: Let S q =o be the total spin 
of the uncharged single-molecule magnet. 

1. If the total spins S q =\ and S q =-i in the ground state 
of the positively (q = — 1) and negatively (q = 1) 
charged single-molecule magnet are equal to S q =±i = 
S q =o — 1/2, then the Kondo Hamiltonian exhibits an an- 
tiferromagnetic exchange coupling, which corresponds 
to the Kondo problem for spin 1/2 impurities or spin 1/2 
quantum dots. 

2. If S q =±i = S q =o + 1/2, then the Kondo Hamiltonian 
exhibits a ferromagnetic exchange coupling that leads 
to a vanishing renormalized transverse exchange cou- 
pling (Ising interaction), in which case the Kondo effect 
is absent. 18 

This point is at the very heart of the Kondo problem. The 
Kondo effect depends crucially on the selection rules for the 
addition and subtraction of a spin on the molecule. This re- 
sult is in contrast to the Kondo effect seen in lateral quantum 
dots where the exchange coupling is always antiferromagnetic 
due to fact that spin states are degenerate in the absence of 
anisotropics 

In the following, we provide a complete and detailed de- 
scription of the Kondo effect in SMMs. Starting from a micro- 
scopic model (Sec. JI]), we derive the effective Kondo Hamil- 
tonian for a SMM attached to metallic leads through tunneling 
barriers (Sec. (TTTI). In Sec. [IVJ we derive expressions for the 
conductance through a SMM for both zero and finite bias as a 
function of a transverse magnetic field and use the new SMM 
N14 to comment the experimental significance of our theoret- 
ical results. Our conclusions are summarized in Sec. IVH 



where U denotes the charging energy, q is the number of ex- 
cess electrons (the charge state of the molecule), and V 9 is the 
electric potential due to an external gate voltage. 20 In the pres- 
ence of an external magnetic field, the spin Hamiltonian of the 
SMM reads 



(q) 

spin 



~A q S qz 



B 



4,9 



(h* ± S q ,+ + h±S q ,-) + h\\S, 



q,z 1 



(3) 



where the easy axis is taken along the z direction and S q ,± = 
S QyX ± iSq yV . The magnetic field components were rescaled 
to h± — gfiB(H x + iHy) and h\\ = gfisHz for the transver- 
sal and longitudinal parts, respectively, where g denotes the 
electron gyromagnetic factor. Note that the transverse mag- 
netic field lies in the xy plane. In this Hamiltonian, the domi- 
nant longitudinal anisotropy term creates a ladder structure in 
the molecule spectrum where the | ± m q ) eigenstates of S q , z 
are degenerate. The weak transverse anisotropy terms couple 
these states. The total spin as well as the coupling parameters 
depend on the charging state of the molecule. For example, 
it is known that Mni2 changes its easy-axis anisotropy con- 
stant (and its total spin) from Aq = 56 /xeV (S q =o = 10) 
to A_i = 43 /ieV (S q =i = 19/2) and A_ 2 = 32 /ieV 
(S q =2 = 10) when singly and doubly charged, respectively. 21 
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FIG. 1: (Color online) Schematic illustration of a single-molecule 
field-effect transistor formed by a single-molecule magnet attached 
to two metallic leads (source and drain) and controlled by a back gate 
voltage. 

The orbital contribution to the SMM energy is given by 



II. MICROSCOPIC HAMILTONIAN 
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The total Anderson-impurity-like Hamiltonian of a system 
formed by a single-molecule magnet (SMM) attached to two 
metallic leads can be separated into three terms (see Fig. [U, 
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The first term on the right-hand side of Eq. (Q]) denotes the 
SMM part, which can be broken into spin, orbital, charging, 
and gate contributions, 



%MM = + ^orbital + U — q eV g , (2) 



where ^ a (Vw) creates (annihilates) electrons in the molec- 
ular orbital state n with spin orientation a and energy e na . 
Here we neglect any diamagnetic response to external mag- 
netic fields. 

The second and third terms on the right-hand side of Eq. 
0]) read, respectively, 
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and 



HSMM-Iead = ~ ^ [^IW^aVW + H.C. 
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and 



where s a = T, k ,k>,*,*< T, a =R,L ^L,a(<w/2)V>fc'. . 

is the lead-molecule tunneling amplitude. The operator 

^la a (^/ecr,a) creates (annihilates) electronic states in the a 
lead (a = L) with linear momentum k, spin orientation a, 
and energy . 



gap of order U between consecutive diagonal elements in Eq. 
©. Since U ^> \t n ,k \, we can assume that the off-diagonal el- 
ements V and W are small perturbations and use a Schrieffer- 
Wolff transformation to decouple distinct charge sectors up to 
terms of order O ( | t n , k | 2 /U 2 ) & 

For simplicity, let us consider the sectors q = — 1, 0, +1 
only, where q denotes the number of excess electrons, and 
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where 



H = A + B, 



(8) 



III. THE KONDO HAMILTONIAN 

The weak tunneling amplitudes between the leads and the 
molecule and the large charging energy cause an effective 
exchange interaction between electrons in the leads and the 
spin of the molecule. Although this interaction at first glance 
seems to have the familiar Kondo s-d form, it is actually un- 
usual because the transversal coupling involves only a sub- 
space of the spin of the molecule. Below, we use perturba- 
tion theory to carefully derive an expression for the effective 
Kondo Hamiltonian of a SMM. 

We begin by considering only the SMM and SMM-lead 
terms in Eq. (Q]). We divide the Hilbert space into subspaces 
corresponding to distinct charge sectors of the SMM. Using a 
block matrix representation, we have 
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where H q is the SMM Hamiltonian for the charge sector q [see 
Eq. (0], while V and W represent the lead-SMM tunneling 
Hamiltonian, Eq. ©. For certain values of the gate voltage, 
the charging energy is compensated and Coulomb blockade is 
lifted. Away from these resonant points, there is an energy 



H-i \ / V 
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Using the similarity transformation H. = e T 7i e _r , where T 
is anti-Hermitian, we have 



n = n + [T,n] + ^[T, [T,n]} 



(10) 



We want to determine T such that B + [T, A] = 0. For that 
purpose, it is sufficient to assume that T has the form 
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with C and V satisfying 

H- 1 C-CH = V 

and 

h v - vn +1 = w, 

respectively. Thus, 

1, 



H = A+-[T,B\ + ()(&), 



(12) 
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where 



cv f + vet o cw-vv 

[T, B) = \ PW f + WD* - Ct V - V f C 
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Note that the neutral sector (q = 0) has been decoupled from 
the charged sectors at the expense of adding two contributions 
of order 0(V, W) 2 to Ho. To specify the form of these con- 
tributions, we use the eigenbasis of H q , namely, 

q (a\H q \(3) q = [H q ] a(3 = SajE^. (16) 
Equations (TT2t and (TT31) can be solved in this representation to 
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In the Coulomb blockade valley, when V g = in Eq. ©, 
the eigenstates of the Hamiltonian Ti q are expressed in terms 
of symmetric and antisymmetric combinations of the eigen- 
states of the S q ^ z operator, namely, 



|s, a) m , = -y= [\m q ) ± | - m q ) ] , (21) 



q=-1 



q=0 



q=1 



FIG. 2: (Color online) Energy spectra of consecutive charging sec- 



tors of the SMM Hamiltonian (q ■■ 



-1,0,1). A q 



(<z) 

I a,,m q 



, q is the 
denote 



tunnel splitting due to the in-plane anisotropy and E t 
energy eigenstates corresponding to symmetric and antisymmetric 
combinations of the ground eigenstates | ± s g ) of the longitudinal 
component of the SMM total spin S qiZ . 



and 



o(<x\V\, 



7?(°) 
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respectively, where [V] a p = _i(a|V|/3) and [W] a p = 
o(ce| W|/3)_i . This allows us to write the following matrix 
elements for the neutral sector: 

= (.9, 



E. 



and 



oHvt C |/J) = E^%^ 



7 ^7 



E 



(o)- 



(20) 



0, 1, 2, . . . , S q if So is integer and m q 



with ra g - w, x,^, ,x ^ 

~, ~, . . . , S q if S g is half-integer. 23 Note that since the ma- 
trices W and V represent the addition (subtraction) of an 
electron to (from) the SMM, only q = states that differ 
in spin projection by 1 are coupled through Eqs. (IT9b and 
(l20l) . Conservation of angular momentum upon electron tun- 
neling requires that the intermediate states in the q = 1 or 
q = — 1 sectors obey |5±i — Sb| = 1/2. For instance, 
when the SMM total spin is lowered by the addition or sub- 
traction of an electron, longitudinal spin components satisfy 
ra±i = S±i = S - 1/2 and -m ±1 -S ±1 = -S + 1/2. 
Therefore, we define the spin states of the ground state to be 
= \m q ) and \[ q ) = \—m q ). The corresponding eigenen- 
ergies are represented in Fig. [5] 



Let us first consider intermediate states involving the q = 1 
sector. This situation corresponds to a small positive gate volt- 
age (V g > 0). Equation (|T9l) can be used to write the ma- 
trix elements of the reduced Hamiltonian of the q = sec- 
tor in terms of a product of energy denominators and matrix 
elements of W. Using the {| | [) q } basis and insert- 
ing complete eigenvector sets for each sector, i.e. 1^ = 
\ s )gq( s \ + \a) qq (a\,we find 



(T \VW^+Wtf\ |)o 
o(l \VW* + WDi\ I)o 



and 



where 



and 



(T \s) 00 (s\V\s) 11 (s\W^\ T)o + o(T |a) 00 (a|©|5)ii( 5 |W t | T)o 
+ o(T \s) 00 (s\V\a) 11 (a\W^\ Do + o(T |a) 00 (a|^|a) 11 (a|W t | T)o + H.c, 

+ oU |5>oo(5|P|a)ii(a|Wt| Do + oU |a)oo(a|»|a>ii(a|Wt| Do + H.c, 



(T IpwUw^I Do = o(T Do + (o(l 1^1 Do) 1 , 

o(T \VW*\ Do = o(T \s)oo(s\V\s) 11 (s\W^\ Do+ o(T |a>oo<a|I>|a>ii<a|W + | Do 

+ o(T |5) 00 (3|P|a)u(a|W t | |)o + o(T |a>oo<a|2^|s>n<s|W t | Do, 

o(I Do = o(l l^oo^l^lDn^l^l Do+ oU I^ooHPI^ii^IW+I Do 

+ oU |5)oo(5|»|a)ii(a|W t | t>o+ oU |a>oo<a|^|a>n<a| 1 Do- 
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Using Eq. d and the fact that E { + 1] = E { + 1] + A lt (T |*>o = oU l s )o = o(T |o)o = V>/2, and (| |o) = -1/V2, we 
find, after some algebra, 



and 
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where A s = E^ - E^, A sa = - E^\ A a 



(27) 



(28) 



(29) 



(30) 



£i 0) - E { a +1) , and A as = £i 0) - E { s +1) . To calculate the charge state q = 1. Similarly, we define 
matrix elements in Eqs. d27ll . (f23b - (t24b . and d25T> we use the 
following definition for the operator W that originates from 
the lead-SMM Hamiltonian in Eq. ©: 



with n being an occupied molecular orbital for a SMM in the 

(32) 



a,/c,cr,n 



o,/c,<T,n 
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A. First set of spin selection rules 

We will now consider the case when adding or subtract- 
ing an electron always decreases the SMM total spin, namely, 
S q =±i = S q =o — 1/2. This selection rule can be enforced 
through the adoption of the following matrix elements: and 
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l(T IV'Ll T)o = Sm , mi -<T $a,i, (37) 



0\T |Vw| T)l -^m 0) mi-(T^,i! 
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id \i>t,\i)o = s. 



mo, — mi—a ^cr,| 



(38) 



(For the sake of simplicity, we will assume that mo, mi > - 
and a = ±| hereafter.) Then, substituting Eq. (l33t and Eq. 



(f34l) into Eqs. ([27]), (|23]), d24j) and (ESJ), we arrive at 
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Using a similar procedure, we also find 

(T |C+V + V + C| T) = E E 4£<]a^^U,a + E E jta7Z'(^'i^ka+^,a'4ua), ^ 



/c,a,n k' ,a' ,n' 



/c,a,n k' ,a' ,n' 
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where 
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and 
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17 + A 



/7- A 
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A±i A±i 



(49) 



(50) 
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I 



In order to arrive at each matrix element expressed in Eqs. account the leading terms in the small ratio A q /U. For in- 
to (l42li and Eqs. d43t to d46k we have only taken into 
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stance, in Eqs. (l39t and (l4Qb we have neglected spin-flipping 
terms which carry an amplitude smaller than the direct, non 
spin-flipping terms by Aq/U. Note that the longitudinal ex- 
change coupling is positive, i.e. antiferromagnetic. Moreover, 
since U ^> Ao, A±i, we find that the exchange couplings are 
strongly anisotropic, with 

|j-L(±l) I w ^l/(±l) (51) 
\ J k,a;k',a'\ ~ jj J k,a;k',a" v Ji ^ 

Finally, introducing pseudospin operators that act solely on 
the q = sector of the SMM, namely, 

= I T)oo(T I - I l)oo(l I, (52) 
= | T>ooU I, (53) 
Sq = IIWII, (54) 



and assuming, for brevity, that Ai = A_i, we can write the 
effective Hamiltonian of the q — charge sector as 



Ho = Ho + [ J k,a;k>,a>^0 (^a^a' " ^, a M,a') ~ Jk,a;k>,a> ( S M|,<A'T,a' + ^0 ^t^Jk'^a' 

k' ,a' k,a 

-3k,a;k',a' (^,a^'T,«' + ^U.a^k'Ua' 



(55) 



where the last term in Eq. (f55b is a scattering term that does not affect the dynamics of the SMM and can neglected. The effective 
exchange coupling constants that appear in Eq. d5Sb are given by the following expressions: 



(56) 
(57) 
(58) 
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The longitudinal exchange coupling is positive, which indi- 
cates an antiferromagnetic Kondo exchange. We have ne- 
glected the dependence of the hopping matrix elements on the 
SMM molecular orbital number n. This is justified when the 
addition or subtraction of an electron brings the molecule to 
the electronic ground state of the particular charge sector. In 
this sense, only one orbital state can be filled (emptied) when 
an electron is added (removed). 

A diagrammatic representation of the longitudinal and 
transverse exchanging interactions is shown in Fig. [3] These 
diagrams differ from the usual Kondo effect in the sense that 
the transverse, spin-flipping interaction in a SMM requires a 
quantum tunneling of the total magnetization during its virtual 
state. 



B. Second set of spin selection rules 

If we consider the case where adding or subtracting an 
electron always increases the total spin in the SMM, namely, 
S q =±i = S q =o + 1/2, we have to modify Eqs. ([35]), ([36]), (l37lh 



and (1351) by adopting instead the following matrix elements: 

o(T IVwl T)i = S m0imi+a S a ^, (59) 



and 



oil \tpmr\ l)l 



-mo,— rai+cr^crl 5 



l(T |^LlT>o = * i 



'rao,mi+cr0<Tt 5 



iU lvLU>o = *- 



(60) 



(61) 



(62) 



Using these selection rules, we arrive at a Kondo Hamilto- 
nian with exactly the same form as that in Eq. (l55l) . The ex- 
pressions for the exchange coupling constants are the same as 
those given in Eqs. ([56]) , d57b , and d58b , except for the longi- 
tudinal coupling J% a . k , a ,, which changes its overall sign and 
becomes negative, signaling a ferromagnetic Kondo exchange 
interaction. However, the strong anisotropy remains, with the 
bare longitudinal coupling dominant over the transversal one. 
We can summarize our results so far by stating that 
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S = m 




(AS z =0, As z =0) 

U 



(a) 



(AS Z = -2m ,As z = +1) 
„ t 2 



u 



2*1 



(b) 



FIG. 3: (Color online) Diagrams representing the (a) longitudinal 
and (b) transversal exchange interactions between the SMM magne- 
tization and the electrons in the leads. 



1. HS q=±1 =S t 



q=0- 



1/2 then J^ a;fe>/ > Jj£ fc , jfl , > 



(antiferromagnetic exchange coupling); 



2. if S q =±i = S q 



- 1/2 then J } 



k,a;k' ,a' 



(ferromagnetic exchange coupling). 



J k,a;k',a' ^ 



IV. CONDUCTANCE AND THE KONDO EFFECT IN A 
SMM 

In order to evaluate the conductance of the SMM subjected 
to the Hamiltonian of Eq. (l55k we make use of the stan- 
dard poor man's scaling to renormalize the effective exchange 
coupling constants J z and J x and the ^-factor. We start by 
calculating the renormalization flow at the points where the 
Kondo effect is observable at zero bias, namely, where the 

tunnel splitting vanishes: Ao (Vi^) = 0- For half-integer 

spins, it is reasonable to assume that ^> Tk, except for 
the first zero (1 = 0). The total Hamiltonian reads 



e m + -?y(ftlE+ + /i±E-) 



Wtot = 

m 

+ ^2£,k1pks1pks +^0 



(63) 



where e m is the eigenvalue of \m) for h± = and, due the 
Knight shift, r\ = 1 — poJ ± /2, with po denoting the density of 
states of the itinerant electrons at the Fermi energy. 19 We do 
not include the Zeeman term for the itinerant electrons in Eq. 
(l63l) because at finite values of one has to cut the edges 
of the spin-up and spin-down bands in the leads to make them 
symmetric with respect to the Fermi energy. 19 We call D the 
resulting band width. 

The Hamiltonian of Eq. (l63l) remains invariant under renor- 
malization group transformations (see the appendix for de- 
tails). Using Eq. (??), we obtain the flow equations 



dJ 1 



= 2p J^J Z 



(64) 



and 



dJ z „ / _ i \ 2 



^_ .2 j± JZ 



(65) 



(66) 



where ( = \n(D / D) and D is the rescaled band width. 
Dividing Eq. (l64l) by (f65t and integrating by parts gives 

(J z ) 2 — (J^) 2 — C 2 , where C is a positive constant. 24 

We have to distinguish between two cases: J z is either pos- 
itive or negative. If J z is positive, the exchange coupling 
constants remain antiferromagnetic during the flow but the 
exchange interaction become increasingly isotropic. Solving 
Eqs. (f64|) and (f65l) yields 



1 



2poC 



arctanh 



= In 




(67) 



The solution for J x is determined by J z = y ( J^) 2 + C 2 . 

The flow of r] is shown in Fig. [5j The flow stops at D w 
ou ~ T > Tk- In the antiferromagnetic case the Berry-phase 
oscillations get strongly renormalized by the scaling of the 
Knight shift. 

Since \J Z \ ^> J ± , when J z is negative, the transverse ex- 
change coupling J 1 - renormalizes to zero. Therefore, in this 
case the Kondo resonance cannot form and the interaction be- 
comes Ising-like. The interesting feature of J 1 - = is that the 
Knight shift vanishes. 



A. Linear conductance 

In order to calculate the linear conductance through the 
SMM we use the following well-known expression for the 
weak coupling regime (Tk <C T):— 



G(T) = Go 



duj 



duo 



16 



\A(uo) 



(68) 



where Go is the classical (incoherent) conductance of the 
molecule, df / duo is the derivative of the Fermi function, and 
A(uo) is the transition amplitude. At the end of the flow, the 
transition amplitude can be calculated in first-order perturba- 
tion theory as 



c 



(lo/T k ) 2 "« c 



{u)ITkY PoG - 1 



(69) 



The Knight shift is related to the transition amplitude by 

P0 A DW 



1 



(70) 



By making the substitution uj — > T into Eq. (1691 ). one finds 
that the linear conductance diverges when T — > Tk, signaling 
the onset of the Kondo effect. Since C > 0, the singularity 
in Eq. (l69l differs from the usual logarithmic behavior found 
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for isotropic exchange interactions. We note, however, that in 
reality, the conductance does not diverge but is rather strongly 
enhanced near the Kondo temperature. The poor man's scal- 
ing breaks down near the Kondo temperature and more ac- 
curate nonperturbative methods, such as the density matrix 
renormalization group, 25 have to be employed for obtaining 
a quantitative description of the conductance dependence on 
temperature. 

Using Eq. d68k we get for the linear conductance 



G(T) 
Go 



16 



(71) 



which has the same functional form as the result of Ref . |2g 
for the resistivity of bulk metals in the presence of Kondo im- 
purities. 

All the zero points of the Berry phase oscillation are 



rescaled by the ^-factor renormalization: = h± } /rjp. 
Thus, the zero points become dependent on the contributing 
states \m) and | — m). This result indicates that the period of 
the Berry phase oscillations becomes temperature dependent 
at T > Tk (see Fig. [4]). This fact allows us to conclude that 
the scaling equations can be checked experimentally by mea- 
suring the renormalized zero points of the Berry phase. Fur- 
thermore, due to the scale invariance of the Kondo effect, the 
distance between the zeros should follow a universal function 
of T/T K (see Fig.g). 



,(0 



B. Nonlinear conductance 

Let us now study the conductance for nonzero bias, 7/0. 
For this purpose, we first express the current flowing through 
the SMM in terms of tunneling rates and lead occupation fac- 
tors, namely, 



dE 



r^r 



Ra 



La 



Ra 



p(E) [f L (E)-f R (E)}, 



(72) 

where p(E) is the energy dependent density of states, Ylo- 
(TR a ) is the escape rate for the left (right) lead, and jp (/r) 
is the Fermi function for the left (right) lead. For the sake 
of simplicity, we assume TLa = ^Ra — T. Since, at low 
temperatures, 



f L (E)-f R (E) 



if Ep - eV/2 <E <E F + eV/2, 
otherwise, 

(73) 

where Ep is the Fermi energy of the leads, we get the follow- 
ing expression for the differential conductance: 



- S(E -E F + eV/2)) . (74) 



Consider the situation where one moves from zero point 
to the magnetic field value b± = + Abj_, where 
Ab±_ = Ah±/rj. If \eV\ <C A (6_l) < T K , the transmission 




Hy[T] 

(b) 




D 



0123456 
Hy[T] 




D 



2 3 4 5 6 
Hy[T] 



(C) 



(d) 





Hy[T] 
(e) 



Hy[Tl 
(f) 



FIG. 4: The graph shows the temperature dependence of the zeros for 
the Berry phase oscillation as a function of the transverse magnetic 
field for the tunnel splittings between |4) and | — 4) states (So = 4) 
and the following temperature values: (a) T /Tk = 1.5, (b) T /Tk = 
1.6, (c) T/T K = 1.7, (d) T/T K = 1.8, (e) T/T K = 1.9, and (f) 
T/T K 2.0. 



amplitude is well approximated by Eq. (1691) . On the other 
hand, for |eV| ^> Tk ^> Ao(6_l), the transmission amplitude 
is given by A e y = J^y. For the case \eV\ ~ Ao(&_l) ^> Tk, 
we can expand ^£>~ max { T a (6^)} u P t0 secon d order in per- 
turbation theory at the end of the flow, yielding 



4 



L 



D-eV/2 T± jz 



D+eV/2 



J± +p J~ ) J~ ) ln 



uo - 
- D 



eV/2 



D + eV/2 



(75) 



where the integration limits account for the asymmetric cut of 
the bands (see Fig. [7]). For \eV\ ~ A (b±) > T, the renor- 
malization flow stops at D = Ao(bj_). Substituting Eq. (1751) 
into Eq. ([74b and setting Ep = 0, we obtain the differen- 
tial conductance up to third order in J ^ for both positive and 
negative biases eV, 



G_ 

Go 



T ^Aj± 2 J z ~ In 



A 



eV 



(76) 



10 




-0.5 
-1 







10 



T/T v 



FIG. 5: The renormalization of the g-factor due to the Knight shift. 
As an estimate, we use poJ± — poJ z = 0.15 where po = 9.45 x 
10 20 J" 1 (seeRef.H3). 



which agrees with the corresponding expression found in Ref . 
|26[ Equation ([76t represents the conductance of the SMM as 
a function of temperature and tunnel splittings between the 
|4) and | — 4) states. In Fig. [6] one can see how the conduc- 
tance depends on the temperature and takes a minimum value 
for those points where the tunnel splitting is zero, which cor- 
respond to the zeros in Fig. [4] Note that our conductance 
formula depends only on the bias voltage and not on the indi- 
vidual chemical potentials of the leads, i.e. our conductance 
formula is gauge-invariant, which is a result of the asymmetric 
band cutting shown in Fig. [7] 

The two split Kondo peaks appear at \eV\ = Ao(6_l). Thus 
the distance between the two peaks oscillates with the mag- 
netic field, following the renormalized periodic oscillations of 
the tunnel splitting Aq(6_l). 



V. CONDUCTANCE AND THE KONDO EFFECT IN A 
SMM IN THE COTUNNELING REGIME 




FIG. 6: Plots showing the temperature dependence of the conduc- 
tance as a function of the transverse magnetic field for the tunnel 
splittings between states |4) and | — 4) for the following temperature 
values: (a) T/T K = 1.5, (b) T/T K = 1.6, (c) T/T K = 1.7, (d) 
T/T K = 1.8, (e) T/T K = 1.9, and (f) T/T K = 2.0. We used a 
bias voltage of V b = V L -V R = 2.7 pV . 



At zero temperature, the current in a single electron tran- 
sistor can be understood as a sequential process of single 
electrons tunneling in and out of the left and right electrodes 
through the SMM, where the transport channel is in between 
the electrochemical potentials of the source and drain reser- 
voir. The electron transport through the SMM is also possible 
for any off-resonant energy, which is commonly called the co- 
tunneling regime, but it is expected to be very small compared 
with the sequential tunneling. The cotunneling contribution 
can be calculated by Fermi's golden rule in second-order per- 
turbation theory, i.e. 



Wi^f = — 



2tt 



(f\H\m)(m\H\i) 



Em, — E 



f 



ill) 

and plays a dominant role whenever the sequential tunneling 
is suppressed, i.e. (f\H\i) = 0. In contrast to the Kondo 
effect, where the tunnel coupling t to the leads is very large, 
in the cotunneling regime t is typically small and thus it is 
necessary to apply a gate voltage V g in order to get close to the 



left 
lead 



right 
lead 



lb) 



Jeft 
lead 



left 
lead 



right + q_ 
lead | 



■eV/2 
-eV/2 



FIG. 7: Diagrams showing the asymmetric cut of the left- and right- 
contact bands when a finite bias is applied. 



resonance condition of for example the q = 1 charged state. 
Thus we need to take only the electron scattering into account, 
thereby neglecting the hole scattering contribution. Using the 
incoherent spin states for temperatures around 1 K, 17 we can 
calculate the cotunneling contribution by means of Eq. (1771) 
in the following form 



2tt 

h 



oU \n tot \s)(s\H tot \ T) 



U/2 - Ai/2 



o(l |W t ot|a)(a|Wtot| ^° 



Z7/2 + Ai/2 



p(E), (78) 
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where U = U + V g . Using Eq. (ED in Eq. CU) we get 



2tt 



o(l \n\ i)oo(T |w| T)o 



17- Ai 



oil \H\ |)oo(T |W| T)o 



(79) 



After some simple algebra, we arrive at the following expres- 
sion for the tunnel rate process, 



2tt 



2Ai 



U 2 -Al 



(80) 



where t is the lead-molecule tunneling amplitude. Comparing 
Eq. (f80b with JjT a -k' a' f rom Eq. (l58l) we see that the second 
contribution to Fermi's golden rule is a function of the cou- 
pling constant, which is to be expected since the Schrieffer- 
Wolff transformation is a perturbative approach of the second 
order. The advantage of using the Schrieffer- Wolff transfor- 
mation is that you can apply the formalism of renormalization 
theory to get a better description of the physical system near 
the Kondo temperature. 

Focusing on the N14 single-molecule magnet, we use 
Eq. (l80l) to calculate the total cotunneling rate between states 
I |) = 1 4) and | j) = | — 4) that will contribute to the current 
flowing through the single-electron transistor, 



4,-4 



2tt_^A| 

h 



/. 



eV/2-U dE 
eV/2-0 (£ 2 "A?)2' 



(81) 



where we integrate over all initial and final states that are 
available within the range of the bias voltage V. Performing 
the integration in Eq. (l8ll yields 



4,-4 



47TI 4 



+ 



U - eV/2 



U + eV/2 



(U - eV/2) 2 - A\ (U + eV/2) 2 - A 2 
(Ax + eV/2) 2 - U 2 



2Aj 



In 



(Ai - eV/2) 2 - U 2 



(82) 



The total current flowing through the single-molecule magnet 
can be calculated in terms of the density matrix by using the 
master equation. Following the same procedure as in Ref.[l7[ 
we obtain the coupled differential equations 



P4 



Ai 

h 



274, 



•74,-4 



(p_4 - Pa) 



(83) 



and 



P-4 



Ai 

h 



2 74,-4 



v?/h* 



• 7I-4 



(P4 - P-a) 



+W-^pA — W^-4P- 



(84) 



where 74,-4 is the incoherent tunneling rate from the lead to 
the molecule. Solving the set of differential equations for p 4 



and p_4 for the stationary case, we obtain the current flowing 
through the SMM for the case where the source and the drain 
leads are oppositely spin-polarized, as described in Ref . [TtL 



I = e W^-aP-a 



2 e 74 _4 A? W4-4 



W4,-4(^ 2 + 74 2 -4^)+ 4 74 2 -4 A r 

(85) 

Figure [8] shows the cotunneling current as a function of the 
transverse magnetic field. Interestingly, the current is sup- 
pressed at the zeros of the tunnel splittings Ai and Ao, exactly 
as in the sequential tunneling regime. 17 



Log 10 I 




100 



FIG. 8: Graph showing log 10 / versus the transverse magnetic field 
for V b = V L - Vr = 4 x 10" 3 eV, V g = O.OleV and 74,-4 = 
10 12 s _1 . One can see that the current is suppressed at the zeros of 
the tunnel splittings Ai and Ao. 



VI. CONCLUSIONS 

The main contribution of this paper is to show how the to- 
tal Hamiltonian of a single-molecule magnet transistor can be 
mapped into the Kondo Hamiltonian by means of a Schrieffer- 
Wolff transformation. While the derivation of the effective 
Kondo Hamiltonian in other contexts, such as quantum dots 
and ordinary single (nonmagnetic) molecule coupled to leads 
is well-known, the case is different for SMM. The dominant 
Kondo effect is unusual for a SMM, since it involves a pseu- 
dospin of the molecule rather than its total spin. We show 
that if the total spin of the molecule is reduced (increased) for 
the charged states, the Kondo Hamiltonian exhibits an anti- 
ferromagnetic (ferromagnetic) coupling, leading to screening 
(anti- screening) of the total spin of the single-molecule mag- 
net. In the case of antiferromagnetic coupling the renormal- 
ization leads to a Kondo effect, i.e. the conductance through 
the SMM exhibits a resonance at the Fermi energy. In the 
case of the ferromagnetic coupling the transverse exchange 
is renormalized to zero, in which case the Kondo resonance 
is absent. This result is in contrast to the case of the Kondo 
effect in lateral quantum dots exhibiting only antiferromag- 
netic exchange coupling, which is due to the fact that all the 
spin states are degenerate in the absence of anisotropics. 19 The 
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standard Kondo screening of the molecule magnetization by 
itinerant electrons in the leads is a very weak effect in this 
context given the large spin of a SMM (its onset is therefore 
likely to occur only at exceedingly small temperatures, inac- 
cessible to current experiments). 

A careful derivation of the effective Kondo Hamiltonian 
shows that the strong dependence of this phenomenon on 
the amplitude and orientation of an external magnetic field. 
The strong uniaxial magnetic anisotropy of the SMM com- 
bined with the weaker in-plane anisotropy creates a pseudo 
spin 1/2 involving states with opposite magnetization orienta- 
tion. A transverse magnetic field modulates the tunnel barrier 
between these states through a Berry phase interference ef- 
fect. That, in turn, modulates periodically the Kondo effect in 
SMMs. 

We have calculated the conductance of the single-molecule 
transistor in the presence of the Kondo effect by using the 
standard poor man's scaling approach. We have show that 
in the case of antiferromagnetic Kondo exchange coupling 
by applying a transverse magnetic field to a single-molecule 
magnet with a large full- or half-integer spin S > 1/2 it is 
possible to topologically induce or quench the Kondo effect 
of the conductance of a current through the single-molecule 
magnet that is sufficiently well coupled to metallic leads. We 
have also shown how the zero points of the Berry phase oscil- 
lation become temperature dependent above the Kondo tem- 
perature and how they change direction within the plane (Fig. 



The latter indicates that the parity of the Berry-phase 
oscillations 7,8 changes from an integer spin S = 4 to a half- 
integer spin S = 7/2. We have also shown how this motion 
affects the temperature dependence of the conductance (see 
Fig. [6]). Interestingly, the maximum value of the conductance 
encircles the zeros of the Berry phase oscillation, providing a 
mechanism for establishing the location of these zeros when 
the orientation of the molecule symmetry axis with respect to 
the metallic contacts is not known. We illustrate these fea- 
tures of the conductance of a SMM using as an example the 
new single-molecule magnet N14. In our view, due to its large 
ground state tunnel splitting this is currently the best SMM 
available for the experimental observation of the Berry-phase 
oscillations of the Kondo resonance. 
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